Binding between endohedral Na atoms in Si clathrate I; 
a first principles study 
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We investigate the binding nature of the endohedral sodium atoms with the density functional 
theory methods, presuming that the clathrate I consists of a sheaf of one-dimensional connections of 
Na@Si24 cages interleaved in three perpendicular directions. Each sodium atom loses 30% of the 3s^ 
charge to the frame, forming an ionic bond with the cage atoms; the rest of the electron contributes 
to the covalent bond between the nearest Na atoms. The presumption is proved to be valid; the 
configuration of the two Na atoms in the nearest 8124 cages is more stable by 0.189 eV than that in 
the Si2o and Si24 cages. The energy of the beads of the two distorted Na atoms is more stable by 
0.104 eV than that of the two infinitely separated Na atoms. The covalent bond explains both the 
preferential occupancies in the Si24 cages and the low anisotropic displacement parameters of the 
endohedral atoms in the Si24 cages in the [100] directions of the clathrate I. 

PACS numbers: 61.50.Ah 
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INTRODUCTION 

The understanding of the mechanism of cohesion of 
condensed matter is essential for solid state physics. Sil- 
icon clathrates are compounds with endohedral atoms in 
the cages of the host frame network and expanded phases 
of diamond type silicon crystal. Cros et al fl|. inspired 
by the structure of the clathrate natural gas hydrates, 
have first synthesized silicon clathrate I containing Na 
atoms. Group 14 clathrates I have been successively syn- 
thesized only when alkaline 0, H, 0] or alkaline earth 
metal atoms Q or CI, Br, or I in group 17 atoms @ 
are encapsulated into the clathrate cages. The electro- 
negativity differences between these host and guest atoms 
are smaller than those in the ionic crystals. If the host 
and guest atoms have large differences, then the induced 
electron transfer forms the ionic compounds with simple 
structures like NaCl or CsCl type structures. 

To date we have found few reports on the role of the en- 
dohedral atoms in the cohesion of the group 14 clathrates. 
The electron charge transfers, from the endohedral Na 
atom to the frame silicon atoms, have been predicted in 
clathrates I and a partial transfer in a Ba@Si2o clus- 
ter [§] . In clathrate II, a displacement of the guest atoms 
has been predicted to be 0.17 A from the center of the Si28 
cage and been explained the displacement to be due to a 
combination of the Jahn- Teller and Mott transition [loj . 
Brunei et al [llj. have observed the displacement using 
EXAFS (extended x-ray absorption fine structure) analy- 
sis. The Na atom was displaced away from the Si28 cage- 
center toward the center of a hexagonal ring by 0.9 ±0.02 
A 11 1 . Libotte et al [l3|. have calculated the displace- 



ments of the endohedral Na atoms in the clathrate II and 
found the displacements to be 0.456 A from the ah initio 
calculation and 0.91 A from a tight-binding calculation. 
Tournus et al. have observed the displacements to be 1 A 
in the the Si28 cage of the clathrate II Na2@Si34 and 2 A 
in clathrate II Nag@Si34 [l3|. They also calculated the 
displacements of the Na atoms in the Si28 cage as 0.65 A 
from the supercell calculation of the Na2@Si5oH44 clus- 
ter with the periodic DFT calculation. They proposed a 
possibility of the displacements to be due to the Peierls 
or Jahn- Teller effect. 

Recently one of the authors has reported the displace- 
ments of the Na atoms in the two adjacent Si28 cages 
hydrogenated to terminate the dangling bonds of the Si 
atoms on the surface of the clusters |l4|. Each Na atom 
displaced by 0.63 A away from their centers of the cages 
to form a dimer between the endohedral Na atoms. The 
displacements was attributed to the formation of cova- 
lent bond between the endohedral Na atoms. They also 
found the electron charge transfered from the endohedral 
atoms to the silicon atoms. 

So the following questions arise: What is the bind- 
ing between the endohedral atoms in the cohesion of the 
clathrates? Why do not the host-guest combinations 
crystallize into the simple ionic structures? In the fol- 
lowing we use a first principles analysis to address the 
questions through investigating the guest-guest and the 
host-guest interactions in the clathrate I. 

Fig. [1] shows a schematic drawing of the polyhedral 
structure of the clathrate I. The structure is special in 
that it consists of the bamboo like Si24 cages only; the 
cages (white) are arranged in bamboos, with spacing a 



2 




fer from the Na atoms to the cage atoms. Next, we eval- 
uate the cohesion energy of the chained Na atoms using 
the periodic DFT method. 



FIG. 1: Polyhedron structure of the clathrate I Si46 by ex- 
tending the simple cubic unit cell. Two horizontal white bam- 
boos of the polyhedron are a one-dimensional bamboo like 
connection of the tetrakaidecahedron (Si24) cages in the [100] 
direction. The connections, arranged in three perpendicular 
directions with spacing a of lattice constant, forms the black 
voids of the pentagonal dodecahedron. This structure is the 
clathrate I Si46, free of endohedral atoms, consisting of the 
tetrakaidecahedra only. 



of lattice constant, in one-dimensional horizontal direc- 
tion sharing hexagonal rings as the bamboo joints be- 
tween the adjacent Si24 cages. Weaving the bamboos in 
three dimensions with common pentagonal surfaces forms 
voids shown with black polyhedron regions in fig. [TJ Each 
void is a pentagonal dodecahedron, separated in space, 
located at bcc position with a different orientation. All 
the previous papers have identified the existence of the 
Si2o cages in the clathrate I. However we presume the 
structure to consist of only the Si24 cages; fig. [T] shows 
that Si2o cages are merely accidental voids in the weaved 
bamboos in the three dimensions. The voids just corre- 
spond to a cages in zeolites, although the voids in the 
clathrate 1 are far smaller than the ones in the zeolites. 
The accidental voids are predicted to have a minor role 
in the cohesion of the clathrate 1. Although this view on 
the clathrate I structure has been neglected so far, the 
experimental preferential occupancies of the endohedral 
atoms in the Si24 cages [l5, 16, 17 1 and the experimental 
anisotropic displacement parameters support this bam- 
boo model. 

So presuming the clathrate I as consisting of the bam- 
boo structures in the three perpendicular directions, 
we analyze the bonding nature between the endohedral 
atoms in the clathrate. First, we calculate the relaxed 
geometries of the one-dimensional clusters with different 
numbers of the Si24 cages using real-space DFT method 
and show the binding nature between the guest atoms; 
the dimer formation due to the covalent bonding between 
the adjacent endohedral Na atoms and the charge trans- 



COMPUTATIONAL DETAILS 



We perform the real-space DFT calculation for the 
hydrogenated bamboo structures using the generalized 
gradient approximation of Perdew, Burke and Ernzerhof 
(GGA-PBE) [13]. We use frozen core 15^2522/ approx- 
imation for the Na and the silicon atoms and the atomic 
orbitals with valence 3s ^ orbital for the Na atom and 
the valence 3s^3p^ orbitals for the silicon atom each for 
which we use double atomic functions for each orbital. 
No smearing for occupations is applied to the final geo- 
metrical optimization. Since we regard the bamboo clus- 
ters as representing the essential aspects of the clathrate 
I, we add hydrogen atoms to the three coordinated sil- 
icon atoms on the surface of the bamboo structures, to 
mimic both the electronic density of states (DOS) and 
the bonding configurations in the clathrate I with the 
clusters. The hydration on the surface of clusters mim- 
ics almost the same electronic states as in the crystalline 
clathrates. The calculated displacements 13, M\ in the 
hydrogenated double Si28 cages in the clathrate II co- 
incided with not only the experimentally observed dis- 

but also the calculated 
O in the crystalline 



placements 0.9 A llj| or 1 A ^ ^ 

displacements 0.456 A [13] or 0.91 A 
clathrates II. This hydration has enabled the states of the 
dangling bonds on the surface of the bamboo structure 
to shift lower side in energy as will be shown in fig. ID 
This hydration has realized the same features as in the 
DOS's of the clathrate Ba8@Si46 [11. We use an ADF 
code [23, [m , which uses a linear combination of Slater 
type orbitals. To evaluate the cohesion energy of the 
chain of the two endohedral Na atoms in the clathrate I, 
we use a periodic DFT code PHASE [i^] with the norm 
conserving pseudopotentials for the Na and the Si atoms. 
For the periodic DFT calculations, Brillouin zones are 
sampled at F and X point set. Markov, Shah and Payne 
have shown that this set is an efficient k-point set to re- 
move defect interactions in the periodic cells [2^. The 
numbers of planewaves are kept 13,805 at F point and 
16,184 at X for any lattice constant. It corresponds to 
set the cutoff energy to be 20.0 Ry (272.11 eV) at 11.0 A. 
We use the PBE [l8| exchange and correlation function- 
als for the electron correlations for the periodic DFT cal- 
culations. We use the spin unrestricted calculations for 
both the real-space and the periodic calculations with 
the convergence of interatomic forces reduced below to 
within 9.45 x 10"^ H/A (5.0 x lO'^ H/bohr). 
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FIG. 2: The inter-Na distances of (a) four-caged Na4@Si78H6o 
bamboo cluster and (b) three-caged Na3@Si6oH48 cluster, 
where the triangles are the inter-hexagonal distances in the 
bamboo structures. The lines are for visual guidance. 
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FIG. 3: The inter-Na distances of (a) four-caged Na4@Si78 
bamboo cluster and (b) three-caged NasQSieo cluster. The 
line is for visual guidance. 



RESULTS 

The distances betvifeen the endohedral Na atoms in the 
relaxed four-caged bamboo structure SirsHgo are shown 
in fig. [2lja) as an example of the relaxed structures of 
even number of cage clusters. Although the distances 
between the hexagonal rings are almost constant, the 
inter-Na distances A and C are however shorter than 



the inter-hexagonal distances: the inter-Na distances A 
(4.84 A) and C (4.85 A) at the ends of the bamboo 
structure are shorter than the distance B (5.38 A). The 
short distances are induced by a bonding between the 
Na atoms. The sum of the shorter and the longer dis- 
tances is 10.21 10.23 A which almost equals the ex- 
perimental lattice constant 10.19 ± 0.02 A of clathrate 
I Na8@Si46 
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We show the distances in the three- 
caged Na3@Si6oH48 cluster in fig. Mi>) as an example of 
odd number of the cages. The inter-Na distances, which 
are smaller than the ones between the adjacent hexago- 
nal rings, are almost the same for each endohedral atom. 
A balance of forces exists between the central Na atom 
and the adjacent two Na atoms. So the small inter-Na 
distances A and C in fig. ^a.) are induced by the dimer 
formation between the Na atoms. The formation may 
leads to a Peierls distortion in the bamboo clusters. 

For the Peierls distortion of one-dimensional case with 
a free boundary condition, the inter-atom distances at 
the edge are different from those in a periodic boundary 
condition. Since two neighbor atoms at the edges form an 
edge state in Peierls gap H^, their distances are longer 
than the ones of inner inter-atom bonding since they are 
located at the free boundary edge. The present bam- 
boo structures have the free boundary condition. Thus 
the Si-H bonds at the edges do form longer Si-H bond 
distances. The Na atoms just inside the bonds in the 
four caged structure form dimers with their adjacent in- 
ner Na atoms as shown A or C in fig. [SJa). The same 
situation occurs for the three-caged cluster in fig. El^b) . 
Here both the atoms forming the distances D and those 
forming the distance E try to form dimers. However they 
are balanced in force. Thus the length D is almost equal 
to that E. 

The distances between the endohedral Na atoms in hy- 
drogen free four-caged bamboo structure Sirs are shown 
in fig.ISja). A single dimer exists at the center of the bam- 
boo structure. Since both the Na atom pairs at the edges 
have formed the edge state forming the relaxed longer dis- 
tance, the Na atoms just inside the bond have formed the 
dimer. Fig. [Sib) shows the distances between the endohe- 
dral Na atoms in the hydrogen free three-caged bamboo 
structure Sigo- The inter-Na distances are the same for 
each endohedral atom; a balance of forces exists between 
the central Na atom and the adjacent two Na atoms. 
These fig. [2] and fig. [3] indicate that the Peierls distortion 
exists between the endohedral Na atoms in these bamboo 
structures. 

Fig. m shows the molecular DOS of the double caged 
Na2@Si42H36 cluster. The shape of the earlier density 
of states 2g, 27 1 of the clathrate are similar to this den- 
sity of states. The HOMO state is at -3.869 eV and the 
LUMO is -3.703 eV, where HOMO is the highest occu- 
pied state and LUMO is the lowest unoccupied state. The 
HOMO-LUMO gap is 0.166 eV. The magnitude of the 
LDA gap has been 0.177 eV with the electron correlation 
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FIG. 4: Molecular density of states (DOS) of the double- 
caged Na2@Si42H36 cluster. The 3s state of the isolate Na 
atom splits into the bonding HOMO state and the several 
unoccupied anti-bonding states. 

by Vosko, Wilk, and Nusair No experimental band 
gap energy has been given, since HOMO is located at just 
above the gap. The eigenvalue —2.754 eV for 3s state of 
the isolate Na atom splits into an occupied single bond- 
ing state -3.869 eV (ISAl.g) which is the HOMO state 
and several higher anti-bonding states 16B3.U (LUMO, 
-3.703 eV), 28Al.g (5.328 eV) and 26B3.U (10.02 eV). 
The decrease of the eigenvalue from the 3s at —2.754 eV 
to the HOMO edge at —3.869 eV is due to the formation 
of the bonding state between the endohedral atoms. This 
is just the bonding state formation in hydrogen molecule. 
Since the HOMO 18Al.g state is composed of a "ger- 
ade" function, the corresponding electron state gives even 
function with respective to the center of the molecule. 
There is a large forbidden region from the HOMO down 
to the state 9Al.u at —6.593 eV indicating the cluster 
is an insulator with HOMO-LUMO gap 2.724 eV, if the 
double caged cluster has no endohedral atom. For the 
four caged bamboo structure, the HOMO-LUMO gap has 
been 0.255 eV. This corresponds to a Peierls gap of this 
cluster. 

To examine the bonding electron distribution between 
the endohedral Na atoms in the double caged cluster, we 
show the electron density profile in fig. [5] given by 

ApNa-Na = ( Na Na ) + ( O O ) 

- ( Na o ) - ( o Na ) 

= p(Na2@Si42H36) +p(Si42H36) 

- p(Na@Si42H36)-p(Na@Si42H36), (1) 

where open circles represent the vacancies of the endohe- 
dral Na atoms. The coordinates of the last three terms 
are fixed at those of the first term to obtain the difference 
of the electron charge densities. This expression gives the 



interaction electron density between the Na atoms, since 
the net number of atoms is cancelled. We have used the 
spin polarized calculations for all the the terms; the non- 
spin states have been the lowest for the first two terms 
and the spin states with /x^ = 1 have been the lowest for 
the last two terms. We have evaluated the sum of the up- 
spin density and the down-spin density for each structure 
and substituted them into the above equation and show 
the density distribution in fig. [5lja). This figure shows a 
clear covalent bonding density between the Na-Na bond. 
This is formed by the dimer formation. The density is 
due to the bonding state between each 3s ^ valence elec- 
tron in the two Na atoms; this is just like the covalent 
bond formation between two hydrogen atoms. We show 
in fig. the difference density on the hexagonal ring 
located at the bisector plane between the two plus marks 
in (a). There is the finite covalent charge densities on 
the plane. Neither the total electron density nor the par- 
tial charge density due to the HOMO state in fig. |4] has 
shown this type of the covalent bond charge densities be- 
tween the Na atoms. The densities of the bonding states 
have appeared between the dimers in the even number of 
cages. 

To see the spatial distribution of the electron trans- 
fers around the endohedral atoms, we show in fig. [6] the 
difference electron density profile 

Ap = Popt - ^ Patom, (2) 

where popt is the density of the geometrically optimized 
cluster and Patom is the overlapped density of the iso- 
lated constituent atoms. The blank zone corresponds 
to regions with lower densities than 10~^ e/A^ or with 
negative densities. So the contour lines correspond to 
the increased charge ones comparing with the overlapped 
isolated atom densities. The electrons around the endo- 
hedral Na atoms are depleted to the cage silicon atoms 
except for the nucleus positions of the Na atoms. 

We calculate electron transfers from the endohedral Na 
atom to the frame atoms. There have been several meth- 
ods to calculate the transfer. Among them the MuUiken 
charges have been found to depend on the number of the 
linear combination of atomic orbitals for the basis func- 
tions [2^. Voronoi charges, which have been named as 
Voronoi deformation charge VDC, have been found to 
gives reasonable values for the transfer ^291. The trans- 
ferred 3s^ electron from each endohedral Na atom to the 
frame silicon atoms have been 0.320e for the double caged 
bamboo structure. The ionic states also appeared in the 
triple caged bamboo structure: The electron transfers 
from Na atoms to frame atoms have been 0.343e (middle 
Na atom) and 0.297e (edge Na atoms) for the triple caged 
bamboo cluster showing that the remained 3s^ electron of 
the endohedral Na atoms has formed the covalent bond- 
ing states between the endohedral Na atoms as shown in 

fig. El 



(b) 




FIG. 5; (Colour print) The spin unrestricted difference elec- 
tron charge density profiles ApNa-Na given by eq. ((TJ, where 
the densities are plotted on a logarithmic scale, 10~^ x 10"''^/^'' 
e/A^, A'^ = 0-10. The blue lines are higher densities than the 
purple ones. The two plus marks correspond to the positions 
of the endohedral Na atoms. The blank regions correspond 
to the densities to be negative or less than 10~^ e/K^ . The 
density (a) is shown on the plane that intersects the two endo- 
hedral Na atoms and the midpoint between two Si atoms on 
the hexagonal ring shared by the adjacent two Si24 cages. The 
density (b) on the hexagonal ring between the two adjacent 
Na atoms in the Si24 cages. 



Here we evaluate the cohesion energy of the Na chain 
in the clathrate I. For this purpose we evaluate the en- 
ergy using the periodic DFT method with the same type 



FIG. 6: (Colour print) The difference charge density between 
the converged self-consistent electron and the overlapped iso- 
lated atom density. The densities are plotted on the same 
plane as in fig. [SJa) with a logarithmic scale, 10~^ x 10^^^^° 
e/A^, N = 0-10. The green contours are higher than the pur- 



ple ones. The blank regions are lower area than 10 ^ 
including negative densities. 



equation as eq. ([I]); 



s/A^ 



Na Na 



Na 



Na 



ST(Na2@Si46) + £;T(Si46) 

£;T(Na@Si46) - £^T(Na@Si46), 



(3) 



where Et^s are the total energies of each crystal. The 
net number of each kind of atoms is also cancelled in this 
equation. The last two terms correspond to that each 
Na atom is located at infinitely separated positions in 
the clathrate. Therefore this equation enable us to eval- 
uate the cohesive energy of the Na chain in the clathrate 
I. The equation has been derived from the difference of 
the formation energies of each phase by Sawada et al [30| . 
They proposed this equation to evaluate the binding en- 
ergies between the substitutional solute atom and the 
interstitial solute atom in the bcc iron. They needed to 
calculate each energy in supercells as large as possible. 
For our calculation the use of the unit cell is sufficient, 
since we need to calculate the cohesion energy of the Na 
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chain in the clathrate. We calculate four kinds of the 
equation of states for each clathrate in eq. ((S]). Here we 
have assumed the energies of the last two terms to be 
equivalent due to their symmetry. The equation used is 

E" = £;T(Na2@Si46) + Et{Shq) ~ 2£;T(Na@Si46). (4) 

The energy of the first term have been the lowest for a 
spin-polarized state /is = 0.188 and the other terms for 
the non-spin states. The equilibrium lattice constant of 
this clathrate has been 10.1998 A, the shorter inter-Na 
distance has been 5.0915 A, and the longer one 5.1083 
A, showing the difference by 0.0168 A in the [100] di- 
rection. The difference between these two distances is 
smaller than that in the hydrogenated cluster in fig.[2ja). 
This is because of the infinite chain connection of the Na 
atoms in the clathrate I. The cohesive energy E'^ of the 
chain has been 0.104 eV which is finite and attractive, so 
the chain is more stable than the two infinitely separated 
Na atoms in the clathrate I. 

To evaluate the energy gain of the distortion in the 
crystalline state, we calculate the total energy of the 
clathrate in which the two Na atoms are located at the 
centers of the gravity of the nearest Si24 cages of the first 
term in eq. This energy has been higher by 0.00186 
eV, with the shorter inter-Na distance 5.0978 A, than 
the full relaxed clathrate. The shorter inter-Na distance 
in the full relaxed clathrate is shorter by 0.0062 A than 
the inter-gravity distance. This quantity is a significant 
difference in the accuracy of the DFT calculations. This 
also indicates the existence of the attractive interaction 
between the shorter Na atom pairs. 



DISCUSSION 

The endohedral atoms have interacted with the cage 
atoms through the ionic bond and with the nearest en- 
dohedral atoms through the covalent bond. 

We have assumed that the clathrate I consists of 
the bamboo structures in the three perpendicular direc- 
tions. Here, we examine the validity of this assump- 
tion. We have calculated the total energy of the clathrate 
Na2@Si46 in which one of the two Na atoms is located 
in Si2o cage and the other in Si24 cage. We have already 
calculated the energy of the clathrate Na2@Si46 in which 
two Na atoms are located at the nearest Si24 cages. The 
energy has been given as the first term in eq. ([4]). The 
energy of this clathrate has been more stable by 0.189 eV 
than that of the former clathrate; the binding between 
the two Na atoms in the chain is more stable than the 
two Na atoms in the Si2o and Si24 cages. This is an- 
other evidence of the validity of our presumption for the 
structure of the clathrate I. 

The covalent bond charge has existed in fig. [5] between 
the endohedral Na atoms. The validity of our bamboo 



structure model for the clathrate I is supported by ex- 
perimental evidence of the preferential occupation of the 
Ba atoms in the Si24 cages by Yamanaka tt al [l^ . They 
reported that the Ba atoms occupy 0.985 of the six Si24 
cages and only occupy 0.189 of the two Si2o cages. The 
high occupancy is a proof of the existence of the covalent 
bond between the Ba atoms Si24 cages. No explanation 
has ever been given for the origin of the occupancies. 

The present study has predicted the difference of the 
inter-Na distances to be only 0.0168 A in the [100] di- 
rection. No report has been existed for the experimental 
guest displacement in the clathrate I except for the guest 
displacement parameters 3l|, [S^- This is because the 



displacement is too small to be measured. 

The anisotropy of the atomic displacement parameters 
of the endohedral atoms in clathrate I has been reported 
by Chakoumakos et al 33, 3J], Nolas et al [s^ in which 



much smaller amplitudes in the [100] directions were re- 
ported than in its perpendicular directions. The present 
study explains the anisotropy to be due to the constrain 
of the displacements of the Na atoms in the [100] direc- 
tions induced by the covalent bond: the bond constrains 
the displacements between the nearest Na atoms in the 
directions. No explanation has been given for the origin 
of the anisotropies. 

The covalent bond between the endohedral Na atoms 
prevents the atoms from crystallizing into ordered ionic 
structure like NaCl or CsCl and crystallizes into the 
caged clathrate structures. The bond forms beads of the 
Na atoms in the clathrate I or three dimensional net- 
work of the Na atoms with symmetry in the clathrate 
II. This is because the electro-negativity of the host 14 
group atoms is smaller than that of the halogen atoms 
that crystallize into ionic crystals. The smaller electro- 
negativity differences between the host and guest atoms 
allow the guest Na atoms to form both the covalent bond 
between the guest atoms and the ionic bond through the 
charge transfer to the cages. Thus the clathrates are a 
compromised electronic state between the ionic crystals 
and the covalent crystals. 



CONCLUSIONS 

Presuming that the clathrate I consists of the sheaf 
of one-dimensional connections of Na@Si24 cages inter- 
leaved in the three perpendicular directions, we have in- 
vestigated the binding nature of the endohedral Na atoms 
with both the real-space and the periodic DFT methods. 
Each Na atom has lost 30% of the 3s^ charge to the 
frame. The finite covalent bonding charge due to the 
Peierls distortion has existed between the endohedral Na 
atoms in the caged clusters. The cohesion energy has 
been 0.104 eV for the chain in the [100] directions of 
the clathrate I. The presumption has been proved to be 
valid; the clathrate encapsulating two Na atoms in the 
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[100] direction has been more stable by 0.189 eV than 
the clathrate encapsulating the two atoms in the Si2o 
and Si24 cages. This covalent bond has explained the ex- 
perimental anisotropic displacement parameters and the 
preferential occupancies of the endohedral atoms in the 
Si24 cages of the clathrates I. The difference between the 
Na-Na distances has been 0.0168 A. This small magni- 
tude of the displacement coincides with the absence of 
the experimental reports on the guest displacements in 
the clathrate I. The beads of the endohedral Na atoms in 
the directions are due to the the covalent bond between 
the endohedral atoms accompanying the electron charge 
transfer from the endohedral atoms to the cages. The 
covalent bond has explained both the preferential occu- 
pancies of the endohedral atoms and the low anisotropic 
displacement parameters in the [100] directions in the 
Si24 cages of the clathrate I. The beads are just a precip- 
itated state in the regular solution theory. The smaller 
electro-negativity of group 14 host atoms than the halo- 
gen atoms allows the endohedral Na atoms to prevent 
the atoms crystallizing into the ionic crystals and allows 
the atoms to form the covalent bonds with beads of the 
endohedral atoms. 

Computations were performed in part using SCore sys- 
tems at the Information Science Center in Meiji Univer- 
sity and Altix 3700 BX2 at YITP in Kyoto University. 
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